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Abstract 

We study the splitting functions for the evolution of fragmentation distributions and the coefficient 
functions for single-hadron production in semi-inclusive e + e~ annihilation in massless perturba- 
tive QCD for small values of the momentum fraction and scaling variable x, where their fixed-order 
approximations are completely destabilized by huge double logarithms of the form a"* -1 \n 2n ~ a x. 
Complete analytic all-order expressions in Mellin-A^ space are presented for the resummation of 
these terms at the next-to-next-to-leading logarithmic accuracy. The poles for the first moments, 
related to the evolution of hadron multiplicities, and the small-* instabilities of the next-to-leading 
order splitting and coefficient functions are removed by this resummation, which leads to an os- 
cillatory small-x behaviour and functions that can be used at N = 1 and down to extremely small 
values of x. First steps are presented towards extending these results to the higher accuracy required 
for an all-* combination with the state-of-the-art next-to-next-to-leading order large-x results. 



1 Introduction 



The fragmentation distributions, or parton fragmentation functions, Dp (x : /j 2 ) encode the proba- 
bility of a final-state parton p in a hard scattering process to end up in (or fragment into) a hadron 
h which carries a fraction x of the momentum of the parent (anti-) quark q or gluon g. Like their 
even more important initial-state counterparts, the parton distributions of hadrons f p h (x,jj 2 ), these 
quantities include long-distance information and are thus not calculable in perturbative Quantum 
Chromodynamics (QCD). Their dependence on the fragmentation (final-state mass factorization) 
scale /j, to be chosen of the order of a physical hard scale in the scattering process under consider- 
ation, is however calculable via the renormalization-group evolution equations 

JLpDf&f) = £ [P/(a.(,i 2 ))®^V)](x). (1.1) 

Here <g> denotes the standard Mellin convolution in the first arguments, 

[fi{#)®htf)\{x) = f^fibfifzi-,!?), (1-2) 
which is reduced to a simple product by the transformation to Mellin moments, 

fi(N,t?) = [ dxx N - l fi{x^ 2 ). (1.3) 
Jo 

The final-state ('timelike') splitting functions Pj in Eq. (11.11) admit an expansion in power of the 
renormalized strong coupling constant, here normalized as a s = OL s (/j 2 ) / (47t), 

P/^cCsGu 2 )) = fl s pJ (0) (x)+a s 2 pJ (1) (x) + fl s 3 pJ (2) W+ ... , (1.4) 

where we have, without loss of information, identified the coupling-constant renormalization and 
mass-factorization scales. 

A benchmark process for parton fragmentation is semi-inclusive e + e~ annihilation (SIA), 

e + e - y*,Z,H -»■ h+X, (1.5) 

which is closely related to deep-inelastic scattering (DIS), eh — > e + X, via the exchange of a 
virtual photon y*, Z-boson or Higgs particle. In the former (latter) case the four-momentum q of 
the exchanged boson is timelike, q 2 > (spacelike, q 2 < 0). X stands for all hadronic final states 
allowed by quantum number conservation. The cross section for vector-boson exchange can be 
written as [HI 

^ra^e = |(i+cos 2 0)p^x,e 2 ) + |sin 2 0p^x,e 2 ) + |cos0F/(x,e 2 ), (i.6) 

where is the angle, in the center-of-mass (CM) frame, between the incoming electron beam and 
the hadron h observed with four-momentum p, and the scaling variable is given by x = 2pq/Q 2 
where Q 2 = q 2 in SIA. The dimensionless transverse (T) and longitudinal (L) fragmentation func- 
tions in Eq. (11.6b , and the total fragmentation function Fj l = Ffi + obtained by integrating 
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Eq. (11.6b over 0, have been measured at LEP and earlier e + e~ colliders, see Ref. [2] for a general 
overview. The interference of vector and axial-vector contributions leads to the parity-violating 
cosG term. The corresponding fragmentation function Fa does not receive l/x contributions and 
will not be considered in this article; it should be ignored when Eq. (|1.6b is referred to below. 

Up to corrections suppressed by powers of the CM energy y/s = Q, the fragmentation functions 
can be expressed in terms of the fragmentation distributions and coefficient functions C a ,i , 

ff&G 2 ) = £ [c a>i (Q 2 )®D?(Q 2 )](x), (1-7) 

f=q,g 

which are short-distance quantities and calculable in perturbation theory, 

C a ,i(x, a s ) = (8 flT 5i q + o^, dig) 8(1 - x) + a s c^J (x) + a 2 c^j (x) + ... . (1.8) 

Here we have identified, again without loss of information, the scale /j in Eq. (11.1b with the physical 
scale Q, and we have suppressed electroweak charge factors. Besides the quantities in Eq. (11.6b . 
we have included the fragmentation functions for the exchange of a scalar § coupling (like the 
Higgs-boson in the limit of a heavy top quark and massless flavours [0) directly only to gluons 

via §G^ v Ga V , where G^ v denotes the gluon field strength tensor. 

The splitting functions (11.4b and coefficient functions (|1.8b are known (with a minor caveat in 
the former case which is not relevant in the present context) to order and a 2 , respectively, see 
Refs. [@]-|9l and references therein. These results provide the next-to-next-to-leading order (NNLO) 
approximation of (renormalization-group improved) perturbative QCD except for Fl, where the 
third-order contributions to Eq. (|1.8b would also be needed. All these quantities exhibit, in contrast 
to their initial-state and DIS counterparts [fT0"l - fl2"l . a double-logarithmic enhancement at small x, 
i.e., the contributions at order a" are enhanced by terms of the form jc -1 \n 2n ~ a x (the minimal offset 
a depends on the quantity under consideration) which correspond to poles a"/(N—l) 2n + l ~ a after 
performing the Mellin transformation (|1.3b . These terms spoil the convergence of the expansions 
(11.4b and (11.8b already at x < 10 2 (81191 and, obviously, preclude describing particle numbers 
(multiplicities) given by the first moments, N = 1, of the fragmentation distributions Dp. 

At leading and next-to-leading logarithmic (LL and NLL) accuracy, these issues were ad- 
dressed long ago, see Refs. [fT3] - [T5l , by showing that these small-* contributions can be calculated 
to all orders. For example, the leading logarithms (a = 2) of P^ g can be resummed to yield 

P T g (N,a s ) k L I (A f-l){(l + g^) 1/2 -l} (1.9) 

with Ca = "colours = 3 in QCD, which leads to a contribution proportional to ^ l /a^ at N = 1. More 
recently there has been renewed interest in the all-x and N = I evolution of (parton) fragmentation 
functions, see, e.g., Refs. lfT6T - [T8l . In particular, a new method has been developed in Ref. |fT9l for 
carrying out the small-x resummation up to the third logarithms in the standard MS factorization 
scheme not adopted in Ref. [fl4l (see Ref. [fTTl for a more detailed discussion of this issue). 

This article builds upon and extends the results of Ref. lfT9l which were mostly presented in 
terms of perturbative coefficients to order oc s 16 . While this is sufficient for collider applications 
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down to * ~ 10 4 , it does neither fix the first moments nor cover the vastly wider *-range relevant 
for fragmentation processes induced by ultra-high energy cosmic rays, see, e.g., Ref. EOl . Here, 
we provide analytic resummed small-* expressions at next-to-next-to-leading logarithmic (NNLL) 
accuracy which facilitate an all-* 'NLO + resummed' evolution of the fragmentation functions, 
and derive the third terms in the resulting N = I expansion in powers of • v /o£. These results 
include the case of Fl, which was dealt with only at NLL level in Ref. [fT9ll . Furthermore we use 
the approach of Dokshitzer, Marchesini and Salam (DMS) in Ref. lETI . which relates the evolution 
of flavour non-singlet fragmentation and parton distributions, to provide an alternative derivation 
of the results for the 'non-singlet' part of P£ and to extend its resummation to the fifth logarithms 

and the ^/a^ 5 contributions at N — 1 in a manner outlined already in Ref. 112211 . see also Ref. Ii23l . 

2 Formalism of the resummation 

Before we present our results, we briefly recall the formalism for deriving the resummation, which 
is based on the mass -factorization relations and the structure of the unfactorized expressions 

F a ,k =C a7 i®Zj k for a = T, <|>, L and k = q, g (2.1) 

in dimensional regularization (we use D = 4 — 2e). The functions C a j are given by Taylor series 
in 8, with the Z k terms including k more powers in In* than the 4-dimensional coefficient func- 
tions (11.81) . The transition matrix Z T consist of only negative powers of £ and can be written in 
terms of the splitting functions (11.41) and the expansion coefficients [3,, of the beta function of QCD, 
P(oc s ) = — Po a s — • • • with P = -y Ca — § «y. This dependence can be summarized as 

fl s M e- n :P r ,p ; a^~ n+l : ... + Pf, ; ... ; a^E' 1 : P„ r _ l . (2.2) 

Hence fixed-order knowledge at N m LO (i.e., of the splitting functions to P m and the corresponding 
coefficient functions) fixes the first m+1 coefficients in the £-expansion of F a ^ at all orders in a s . 

The small-* expansions of Fj^ (the corresponding relation for Fi is slightly different) read 

F\ a „ e - n+i = f$x- l \n n+£ - l x+F$x-hri n+£ - 2 x + ... . (2.3) 

If the constants up to ^7 316 known for all n and I, then the splitting functions and coefficient 
functions can be determined at N m LL accuracy at all orders of the strong coupling. As observed 
in Ref. [fT9l . the n th order small-* contributions to Ft.§ are built up from n terms of the form 

(A n . k e- 2n+l + 5,a£- 2 "+ 2 + ...Jx- 1 " 2 * 8 , k=l,...,n. (2.4) 

Since the terms with e~ 2n+1 , ... , £ - " -1 have to cancel in sum (|2.1I) . there are n— 1 relations 
between the LL coefficients A n ^ which lead to the constants £ | in Eq. (12.31 ), n — 2 relations 
between the NLL coefficients B n ^ etc. As discussed above, a N m LO calculation fixes the (non- 
vanishing) coefficients of £~", .. . , £~ n+m at all orders n, adding m+1 more relations between 
the coefficients in Eq. (|2.4I) . Consequently the highest m+1 double logarithms, i.e., the N m LL 
approximation, can be determined order by order from the N m LO results. Finally the resulting 
series, here calculated to order a^ 8 using Form and TForm (|24|. can be employed to find their 
generating functions via over-constrained systems of linear equations. The whole procedure is 
analogous to, if computationally more involved than, the large-* resummation in Ref. [25 J. 
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3 Af-space results: splitting functions P£ 

It turns out that the resummed splitting functions in Mellin space can be expressed in terms of 

S = (l-^) 1 / 2 and L = ln(i(l+S)) = -ln(^(l-5)) (3.1) 

with \ = -SC A a s /N 2 and N = N-l. At NNLL accuracy, i.e., resumming the contributions 
a" x~ l \n 2n ~ a x with a = 2, 3, 4 the (flavour- singlet) splitting functions in Eq. (11.11 ) are given by 

+ ^ a s N | (- 1 1 C\ + 6C A n f - 2QC F n f ) ^ (S - 1 + 2$) + 10Q 2 1 (5 - 1) X 

- (51 Cj - 6C A n f + \2C F n f ) \ (S - 1) + (1 1 C\ + 2C A n f - AC F n f )S~ x L 

+ (5C A 2 -2C A n f + 6C F n f )^(S-l)L 2 + (5lC A 2 -l4C A n f + 36 C F n f ) L } , (3.2) 

P Zg( N ) = § P qqW - l^a s N{C A 2 + C A n f -2C F n f ){^(S-l)(L+l) + l] , (3.3) 
P* g (N) = \N(S-1) - ^a s (UCl + 2C A n f -4C F n f )(S- 1 -l) - P^N) 

+ jj^j a s N | ([1 193 - 576 &]C A - 140 C\n f + AClnj - 56C^C F n f + \6C A C F nj 
-A^Cjnj) (5-1) + ( [830 - 576 t&]C} + 96Cln f -SC^n} -20&C A 2 C F n f 
+ 64C A C F nf -96Cjn])(S- 1 -1) + (\\Cl + 2C A n f - AC F n f f(S^ - 1)} , (3.4) 

P U N ) = ^Pg g (N) - ^a s (C A 2 + C A n f -2C F n f )±(S-l+2t > ) 

+ 3£ §7 ( U C a + 13c a"/( c a - 2C F ) +2Clnj - 8 (C A - C^C^ 2 ) (1 - S" 1 ) 

- (48 C| - 45 C|C F - 72 ^ 2 Cj (C A - C F ) - 33 Cjny + 2 c\n) + 48 C^C^ 

- 8 Cjnj) | (5 - 1 + 2 £) + (-54C? + 45 C A 3 C> + 72 £ 2 C A 3 (C A - C F ) + 23 Q 3 ^ 
-28C A 2 n / CF-8(CA-2C F )C F n/) |(5-l)x} (3.5) 

with C/7 =4/3 in QCD and ^ 2 = ^ 2 /6. The respective first lines of Eqs. (13T21 - (1331 ) are the LL (for 
Pgq and P g g, already determined in Ref. lfT3l0 and NLL contributions, the rest represents the NNLL 
terms. Of course, no negative powers of C A and no non-(/V— 1) _1 terms remain when these results 
are expanded in powers of a s . After combination with the LO and NLO splitting functions flU - 
with the l/(/V — 1) , 1 < £ < 3 poles removed to avoid double counting - these results provide a 
combined (7Y= 1 finite) all-x 'NLO + resummed' evolution of the MS fragmentation distributions. 
The crucial step towards deriving Eqs. (13.21) - (13.51) is discussed in Appendix A. 

As discussed in Ref. |fT9l , it is possible to also obtain the next (N 3 LL) contributions to P q T q and 
P q g, due to A„ i = B n \ = . . . = in Eq. (12.41 ) for F F . q and F^q. We have been able to find the 
exact all-order expression also for these contributions. However, the results are considerably more 
lengthy than Eqs. (13.21) and (13.31) and hence are deferred to Eqs. (IB. 1 1) and (IB.2I) in Appendix B. 
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The first moments of the combined splitting functions receive contributions from (13.21 ) - ( 13.51) , 

T(m) 

taking the limit N — > 1 for fixed oc s , and the 'truncated' fixed-order results P with 

^qq - ) ^qg ~ 3 "/ > ^gq ~ ~ 3C ^ ' ^gg ~ ~ T Ca ~ 3 "f 

atN = 1. This leads to the 'NLO + resummed' results 

P q T g (N=l) = ^n f a s - i {nCjtn f -2C A nf + 4C F n}) (2C A a s 3 ) 1 / 2 + 0(a 2 ) , 

/g(iV=l) = (2C A a s ) 1 / 2 - ^(11^ + 2^+ 12^) « s 

+ -f^r ( [ 1 1 93 - 576 C2] c£ - 1 40 C\n f + 4 C 2 /^ 2 + 760 Cj(C F n f 
-mC A C F n}+ UACjnj) (2C A a 3 ) 1/2 + 0(a 2 ) , 
/5 gq(^ V=1 ) = % (^=l) + f^a s ) + 0(a s 2 ). (3.7) 

It is interesting to note that the combination P qq — PT~ + P gq — P^L of the resummed first moments 
(13.71 ) vanishes for C A = C F irrespective of the numbers of flavours n*. Instead inserting the QCD 
values for the colour factors, we obtain for n^ = 5 the numerical series 

p^ q (N=l) £ 0.2358a s - 0.6773a s 3/2 + 0.5880a s 2 , 

Pq g (N=l) = 1.0610a s - 1.5240a s 3/2 + 1.8089a s 2 , 

P g T q (N=l) 0.3071 aj 2 - 0.3059a s + 0.2884a s 3/2 , 

P g T g (N= 1) £ 0.6910a s 1/2 - 0.9240a s + 0.6490a s 3/2 , (3.8) 

with benign coefficients in terms of a rather large expansion parameter with v /a^ ~ 0.34 at /j = Mz- 

In order to illustrate this behaviour we have included the a 2 coefficients resulting from Eqs. (IB.ll) 

Yd) 

and (|B.2I) together with the corresponding truncated second-order splitting functions P ■ . 



4 Coefficient functions for F T and F§ 

We now turn to the coefficient functions ( 11.8b for Fj and F§. Their leading and next-to-leading log- 
arithms, a"x~ l \n 2n ~ a x with a = 1, 2, need to be included in a 'NLO + resummed' approximation 
which is applicable at all values of x and finite at /V = 1 . The corresponding /V-space results, again 
derived from the series expansions in Ref. lfT9l . are given by (using F = S~ 1//2 ) 

C r , g (iV) = g(F-l) + - n ^2^/?{(-555C A 2 + 66C A n / -480C F n / )(F- 1 -l) 

+ (868C 2 + 152 - 336C F n f )(F - 1) + (1 1 C 2 -2C A n f ) [6 (F 3 — 1) 
-5(F 7 -1)] -(132C A 2 + 24(C A -2C F )n / )( J F 5 -l)+384C F n / Fi:| , (4.1) 

CV, q (A0 = \ + \-^N{F- l -l-FL) (4.2) 



and 

C fg (iV) = l + ^C T ^N) + T ^N{(2C^-C A n f + 2C F n f )(F- l -l) 

- (ACl + C A n f -2C F n f ){F-\)} , (4.3) 

C ^ N ) = § (PraM-l) • (4-4) 



Here we have included, besides the resummed 1 / (N—l) terms, the zeroth-order terms in Eq. (11.81) . 
The leading logarithmic contributions in Eqs. (14.11) and (14.31) have been obtained before, in a com- 
pletely different manner, in Ref. IfTTl (where, as in Refs. (HE!, the normalization of C fl)g in SIA 
differs by a factors of two from that adopted in the present article and Ref. [0). We have also been 
able to derive closed ,/V-space resummations of the NNLL contributions, a"x~ l In 2,1-3 x, to all four 
quantities. The corresponding more lengthy expressions can be found in Appendix B. 

The first moments of Eqs. (|4.1I) - (|4.4I) vanish except for contributions in Eqs. (|4.1I) and 
(I4.2I) . Note, in particular, that the zeroth-order contribution (II .8b is compensated in C^ g by the 
resummation of the leading logarithms. The first a s -dependent contributions are due to Eqs. (IB.7b 
- (IB. 10b and the truncated NLO coefficient functions cjty (i.e., the corresponding quantities of 
Eq. (|1.8b with the 1 / (N—l) poles removed) and read 

C T>q (N=l) = -§C T)g (iV=l) + ... 

= l + i^(9Ci + 50C2n / -8(C A -2C F )n/)a s +..., (4.5) 

2 

C^, q (N=l) = -§^ )g (iV=l) + ... = -±%(C A -2C F )a s +... . (4.6) 

The next corrections at N = I may be expected at order a 2 which is beyond our present reach. 

All our results above are given in the standard MS factorization scheme. The transformation to 
another scheme S is given by 

C S (N) = C^lZsiN)]- 1 , P S (N) = (Z S (N)P(N) + V^§^)[Z S (N)}- 1 , (4.7) 
where (suppressing the Mellin variable and the scheme index) C and P represent the matrices 

C= ( C J« °jA and P = (f T f \ , with /> LL = f ° ) (4.8) 



in MS. Eq. (14.7b includes, for Z5 = C, the transformation to the 'physical evolution kernels' K, 

(4.9) 



d / Ft \ _ / ^rr £r<|> \ ( Ft 

d\nQ 2 v i> r\% % y V F 

for the system F = (Fr, F^) of fragmentation functions [for N ^ 1, since det C vanishes for the 
results in Eqs. (14.5b and (|4.6b l. cf. Ref. 11261 . At LL accuracy this transformation, as well as that to 
the massive-gluon scheme of Ref. lfT4l . see also Ref. IfTTl , is of the form 

1 ( 1 -''.S.I.I. Q ( ( 'l.l.- ( '.S.I.I. ' I (4 



-5,LL 



i+^LLV 1+C 



LL 
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with c LL = F — 1 [recall F = S -1 / 2 with S given in Eq. (13.11 )1, c s LL = i (5 _1 — 1) for the scheme 
of Ref. |fl4l and c SLL = for the physical kernels. P LL is invariant for this form of Z5. Since we 
do not know the coefficient functions in the massive gluon scheme beyond LL accuracy, we are 
not able to use the scheme transformation to compare our results to the NLL splitting functions 
in Ref. lfT4l . The transformation to the physical kernels in Eq. (14.91 ), on the other hand, does not 
pose problems. It is interesting to note that, disregarding overall prefactors, the results to NNLL 
(for K§t and K§§) and N 3 LL (for Kjj and Kj^) accuracy include only integer coefficients in the 
expansion in powers of a s . Consequently terms with the logarithm L in Eq. (|3.1I) do not enter their 
closed all-order expressions, in contrast to Eqs. (13.21) - (13.51) for the MS splitting functions. 



5 Coefficient functions for F L 

The remaining SIA coefficient functions with x _1 \n k x contributions are those for the longitudinal 
fragmentation function F F . There is no oc^ term in this case, hence the a" results are required for 
the N' 1 *LO approximation and thus for resummation of the n highest logarithms outlined above. 
Since the NLO coefficient function C^g includes terms up to x~ l In 2 *, the NNLL resummation is 
required to extend the dX\-x 'NLO + resummed' approximation to Fl. The coefficient functions for 
Fl have been calculated so far only to order a 2 (5][7]|. It is, however, possible to derive (at least) 
their highest three a 3 * -1 ln^jc contributions using x — > l/x analytic-continuation (RC) relations 
between DIS and SIA physical kernels along the lines of Ref. (9), see also Ref. (27). 

For this purpose we have calculated the evolution kernels K T and K s corresponding to Eq. (14.91 ) 
for the system (F T ,F L ) with F L (N,Q 2 ) = F L (A^,2 2 )/(a s c}|q(A^)) and its 'spacelike' counterpart 

of DIS structure functions (Fj, Ff) for which all NNLO coefficient functions are known ifTOlfTTl 
l28l . At this order these kernels are expected to be related by 



K T ab {x,a,) = K* b (x,a 



(5.1) 



up to terms beyond the presently required logarithmic accuracy. These relations form a check of 

(3) (3) 

the respective two highest logarithms in ' and L which have been derived by resumming the 
NLO results in Ref. [19J, and provide an over-constrained system of four equations for the hitherto 
unknown coefficients of oc s 3 x~ l \n 2 x in CL, g and a s 3 x _1 \nx in C^ q . Eq. (15.11) leads to 

c £) (AT) = 256 C A C F n f AT 4 - 64 C F C A n f - | C F n f - | nfj N~ 3 

- 64C> C A n f - [| - f C2] C F n f - \ nfj N~ 2 + 0(^ 1 ) , (5.2) 

c { ^ g (N) = UQ%ClC F N- 5 -32C F (^-Cl- 10 C A C F - ^ C F n f - \c A n^ AT 4 

+ 32C F ([19-22^ 2 ]C A 2 -f C A C F + ^C A n f -^C F n/)N- 3 + 0(^ 2 ) . (5.3) 

The resulting NNLL resummation for Fi - which does not involve Eq. (15.21) - can be written as 
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+ T728 |96 (5C}-2C A n f + 6C P n f ) ^(S - \)FL 2 - (23 -42 C A n f 

-32C A C F + \2SC F n f ) 9|(S , -l)F~ 1 i: + 8 (53Cjt + 2$C A n f + 36C A C F 
-114 C F n f ) | (S - 1) F£ - 6 (55 C 2 + 6 C A n f - 16C F n f ) ^(S - 1) F 3 L 

- 24 (1 1 C\ - 2C F n f ) |(5 - 1) F 5 £ - 5 (l 1 C A 2 - 20«/) Q (5 - 1) F 1 L 
+ 4(F 5 — 1) j — 16 (71Q 2 + 10C A n / -48C F n / ) i ((5 - l)F~ l + 2$) 

- 12 (577 Ci - 54Ca/» / -96CfCa + U4C F n f ) (F -1 — 1) +4 (755 C 2 

- U6C A n f + 384 0«y) (F - 1) - 4 (209 Cf + 10 O"/ - 4SC F n f ) (F 3 - 1) } ,(5.4) 

+ m % «s { 1 6 (29 Q 2 + O^/ - 3 C F n f ) ( ± (5 - 1 )F~ 1 + 2) + (737 C A 2 + 58 O"/ 
+ 28800- 1920"/) (F- 1) - (l21Cjt + 26C A n f -4SC F n f ) (F 3 - 1) 

- (209 C 2 + 10 C A n f - 48 0«/) (F 5 - 1 ) - 5 ( 1 1 C\ - 2 C A n f ) (F 7 - 1 ) } 

+ 33P776 §r a s N { - 256 ([104 - 4212 £>] C# - 459 C|0 + 533C|«y - 480 ClC F n f 

- \6Clnj -6QC F nj[C A -3C F ]) (^(5— 1)F _1 + 2) - 368640n / (C A 2 + On / 
-2C F n f ) (|(S - 1)FX) + ([221217 + 101376 £ 2 ] - 77760C|O - 3852Cj/z f 
- 9856^0^- 512 [12CAC}n f + C A C F Tif- 3 C|t$ -92C\nf) 9(F~ l - 1) 

- ([100985 - 1575936 £ 2 ] O? -451584CjO + 234740(7^ - 160032C|C F n / 
+ 19353600 2 "/ - 14524C 2 "/ + 8774400"/ - V2M\6Cjnf) (F - 1) 

- 3 ( 1 1 6329 + 55296 £ 2 + 23668 Cjny + 1924 C^nj-l 2672 cjo 
-55116C%C F n f - 12480 C A C F n} + 16128 Cjn 2 ) (F 3 - 1) - 3 (126471 O? 
+ 55296 ^ 2 + 3800^(5 C A + n / ) + 25344C A 3 C F - 50752 Q^O"/ 

- 9344 OO" 2 -92 16 00%+ 176640^/) (F 5 - 1) -3 (71511 C} 

- 182600% - 1348 C 2 " 2 + 10560C|O - 9664C 2 On / + 6016OO"/ 

- 23040%) (F 7 - 1) + (123541 + 55396C% - 26360% -(lOC A C F nj 
+ 209^0^)480 + 1 1520 Cjnj) (F 9 - 1) + (4411 C\ - H2C 3 A n f - 44 Clnj 
-9600^/(110-2^)) 35(F 11 -l) + 385Cf(HO-2n / ) 2 (F 13 -l)| . (5.5) 

Eqs. (15.41 ) and (15.51 ) complete the results required for the 'NLO + resummed' all-x flavour-singlet 
evolution of the fragmentation distributions (11.11 ) and fragmentation functions (11.6b . It is interesting 
to note that, unlike the corresponding splitting functions and kernels for the system (F F , F^), the 
upper-row elements K Fa of the resummed evolution matrix for {Fr-. Fi) are suppressed by two 
powers of N = N— 1 with respect to their lower-row counterparts K Fa . Also these kernels are found 
to include only integer coefficients after expansion in powers of a s = 0C S / (4n) or a s /7L 
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6 DMS relation and 'non-singlet' resummation of P^ g 

Another connection between the final-state ('timelike') and initial-state ('spacelike') evolution has 
been suggested in Ref. I12TTI . The evolution of the flavour non-singlet fragmentation distributions 
and parton distributions, respectively denoted as /j and /_ l5 is written as 

aipM^Q 2 ) = [Pu(a s (Q 2 ))®Mz°Q 2 )](x) (6.1) 

(the convolution cg> has been defined in Eq. (11.21) above), and the modified splitting functions P u 
are postulated to be identical for the timelike and spacelike cases. Eq. (16.11) has been applied and 
verified for the NNLO non- singlet and gluon-gluon [[8J splitting functions. In the latter case its 
applicability is obvious for Quantum Gluodynamics, = 0, but it is found to hold for all terms in 
the limit C F = 0. In /V-space Eq. (IBTTT) leads to BUI 

d lnQ 2MN,Q 2 ) = P„(N)MN,Q 2 ) = P u (N + od lnQ 2)MN,Q 2 ). (6.2) 

The solution of this equation, obtained by expanding in powers of o, can be cast in the (new) form 

Pa(N) = P u{N)+ i^^L^ [Pu{N)r ) (6.3) 

n=l \ ' 

suitable for all-order considerations. The above relations imply that the difference hP = P\—P \ 
of the time- and spacelike splitting functions at any order is given by lower-order quantities. 

The crucial point in the present context is that the spacelike splitting functions are only single- 
logarithmically enhanced with the (scheme-independent) leading terms H29H30H 

P£g(N) = £f^N- 1 + 2^ (^) 4 jV- 4 + 0(a s 6 jV- 6 ) . (6.4) 

Therefore the above differences and the timelike splitting functions are equal, 8Pgg = P gg , for 
n > £ + 1 at the level of the N^LL small-x double logarithms. The LL, NLL and NNLL contribu- 

T(yi\ 

tions to P gg for C> = are thus fixed by Eq. (16.11) to all orders in a s without any spacelike input 
beyond the NNLO splitting functions [[T2l. and are found to be identical to the results presented 
in Eq. (13.41) . As already discussed in Ref. [|22l . the N 3 LL terms are fixed if Eq. (16.41) is used in 
addition, and finally the N 4 LL corrections require only one additional coefficient, that of a^N^ 3 , 
of Pg g (N). However, this coefficient is presently not available in the literature for the MS scheme. 

The resulting prediction for the N 3 LL contribution to P g T g reads 

P U N )\%ll = T3l24^«s{4([48125-3168C 2 - 15552^]ci-[46 + 2880C2]C A 2 n / 
+ 236C A nf - 8 nj) | (S - 1 + 2$) - ([48473 + 9504 ^ 2 - 25920 £3] Cj - 340 
+ 24 n) - [ 2670 - 4032 £ 2 ] C\n f ) ^ (S - 1 + 2 \ + 2 Z?) + 2 ( [ 20337 + 7920 £ 2 
- 12960 y Cl - [2330 - 2592 £2] C\n f - 156 C A nj + 24nf)^ (S" 1 - 1 - 2\ - 6^ 2 ) 
+ 4([ 1617 - 1584^ 2 ] C\ + [316 - 288 £ 2 ] C\n f - AQC A nj - %nj) ^ (S~ 3 - 1 
-6^-30^ 2 ) + (llC A + 2n / ) 3 ^(5- 5 -l-10^-70^ 2 )} . (6.5) 
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Combined with the first moment of the truncated NLO splitting function, 



-p T W(M n 160 r-2 j/76 ,16 r \ 88 r n ~ 

F gg ( N = 1 ) = ~TT c a + I 27 + T^ 2 J A / _ T F f ' ^ ^ 

this results yields 

'£("=l)£r° = -(¥-Tfc-6fo)ci- (^"f b)Cx^. (6.7) 
Including also the N 4 LL contribution in Eq. (IB.12I) . the numerical result for ny = 5 is given by 

p r (/y=l)| CF=0 0.6910a s 1/2 - 0.5703a s + 0.0267as /2 + 0.4946a s 2 



+ ( 1.0036 + 0.1 5753 fll (3) 



c f=°\ 5/2 , 3 

NLL 



ar + 0(a s J ) , (6.8) 



where 5 gg | nll is the above-mentioned next-to-leading small-* coefficient at order a 4 , cf. Eq. 
(11.41) . in the notation of Eq. (16.41) . i.e., the coefficient of (CaOC s /tc) 4 A r_3 . One may expect it to 
be negative and no larger than a few times the corresponding LL coefficient 2C,^ — 2.4041. An 
explicit determination of this coefficient and its fifth-order counterpart from the results of Ref. 11311 
should definitely be possible, in particular in view of the calculations performed in Refs. [|32l[33l . 
On top of these quantities, and the known NNLO splitting function [|8l and sixth-order coefficient 
in Eq. (16.41) . the aj coefficient in Eq. (16.81) requires the NNLL fourth-order term in Eq. (16.41) . as 
already discussed in Ref. [|22l . which currently appears to be out of reach. 



7 Analytic and numerical results in x-space 

We finally return to the now complete 'NLO + resummed' approximations for the MS splitting 
functions and coefficient functions in Eqs. (11.11) and (11.71) . Their numerical consequences in 
x-space can be obtained either by expanding the above /V-space results to any desired order in 
cc s and using 

k\ 



Jo 



dxx N - 2 \n k x = (-l)S — — r— -j- , (7.1) 

or by subjecting the iV-space expression to a standard numerical Mellin inversion, cf, e.g., Ref. [1341 . 
The former can also be used to identify the closed form of x-space results such as 

X pT(x,a s ) +*P q J(x,O g )| NNLL = |4C A a s + | (HC 4 2 + 2C A n / -4C F n / )a s 2 lni} \j x (z) 

+ { \ (260^ - 23 C A n f ) a s 2 + ^ ( 1 1 Cj + 2C A n f - 4C F n f ) 2 a 3 s In 2 \ } \ J x (z) (7.2) 

+ 9% ( [ 1 34 ~ 72 c l + 23 C\n f - 48 C 2 A C F n f + 4 C A C F nj - 8 Cjnj) a 3 In 2 \ £ J 2 (z) 

and 



xP^(N) - g a.) = -f ^(Cjf + C A i^-2C F i^Klii± J7 2 (z) (7.3) 



32 C f i f~* M ^ „ x„2t„ 1 4 
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with 

z = (32C A a s ) 1 / 2 ln^ (7.4) 

in terms of the Bessel function of the first kind J\{z) and Ji(z), cf. Ref. [|35l . The first line of 
Eq. (17.21) represents the LL and NLL terms, and the next two lines the NNLL contribution. The 
missing ingredients for a completely analytical Mellin (Laplace) inversion of Eqs. (13.21) - (13.51) are 
the inverses of products of powers of the square root S and the logarithm L defined in Eq. (13.11) . 
Both this logarithm, which also corresponds to a Bessel function in x-space (see Eq. (IA.3I) below), 
and the fourth root F in Eq. (14.11) . which corresponds to a hypergeometric function ifTTl . can be 
viewed as 'artifacts' of the MS scheme, as neither remains after transformation to the matrix of 
physical evolution kernels (14.91 ) for the fragmentation functions Ft and F§ and their counterparts 
for the system (Ft, F£) discussed above Eq. (15.11) . 

The result (17.21 ) for Pj„ indicates a general single-logarithmic enhancement of the Bessel- 
function oscillations at extremely small x, where furthermore the contribution of Ji(z) is sup- 
pressed by a factor 2/z with, e.g., z ^ 0.96 In- for a s ~ 0.12 - recall that J\(z) and J2(z) differ 
for large z only by a phase shift of k/2 11351 . The asymptotically dominant a s (a s ln-)^ 2 Ji(z) terms 
arise from the S~ 2(+l contributions to the N^LL terms and have closely related prefactors propor- 
tional to ( 1 1 C\ — ICavi^ + ACfh^ , see Eqs. (13.41) . (16.51) and (IB. Ill) , which seem to point towards 
the possibility of a 'second resummation' . Taking into account the asymptotic z~ 1//2 suppression 
of J n (z), the oscillation amplitudes of these contributions to the LL, NLL and NNLL results behave 
as (ln-^) -3 / 2 , (ln-^) -1 / 2 and (In ~) 1 ^ 2 , respectively. 

The LO + LL and NLO + NNLL approximations, which provide the minimal N=l finite re- 
summations of the respective fixed-order results, are illustrated in Fig. [Qfor the second column and 
in Fig. |2]for the first column of the splitting-function matrix P in Eq. (14.81) . Note that the scales of 
the right and left parts of both figures are related by a factor Ca/Cf = 9/4. The close similarity 
of the curves at small x thus directly demonstrates the approximate 'Casimir scaling' of Pj and 
Pj. The oscillation amplitudes of the former quantities are much larger than those of the latter. 
This behaviour is not due to the LL contributions to Pj which, as for other resummations, are 
numerically small. There is no reason to assume that terms beyond the present NNLL accuracy 
will be small. The Mellin inversion of Eqs. (16.51) and (|B.l II) . however, appears to indicate that the 
amplitudes will stabilize, for most of the wide x-range shown in the figures, after including the 
N 3 LL and N 4 LL contributions. 

Corresponding minimal 'LO + resummed' and 'NLO + resummed' jc-space results are shown 
in Figs. [3] and |4] for the coefficient functions for Ft (left) and Ft (right). In this counting, the 
leading contribution Cf/Cac ll to C^g in Eq. (14.11) is not combined with the 8(1— x) LO term in 
Eq. (|1.8I) . This is consistent with the fact that c LL is a scheme-dependent quantity as discussed 
around Eq. (14.101) . Given the oc s -expanded N 3 LL result in Ref. lfT9l , we expect that the coefficient 
function Ct a , for which Fig. |4] includes the Mellin inverse of Eq. (|4.2I) . will be the first quantity 
determined at an N = 1 finite 'NNLO + resummed' accuracy. The opposite is true for C^g, the least 
stable quantity at the present accuracy as shown in the right part of Fig. |3] Here, as for P?, the 
highest five logarithms need to be resummed for an all-x combination with the NNLO results. 
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Figure 1: The timelike gluon-quark and gluon-gluon splitting functions (multiplied by x), shown 
for a very wide range of the momentum fraction x at a typical value of the strong coupling constant 
a s . The all-x (minimal A^= 1 finite) 'LO + resummed' and 'NLO + resummed' approximations are 
compared to the corresponding LO and NLO results valid only at large x, e.g., x > 10 for NLO. 
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Figure 2: As Fig. \T\ but for the timelike quark-quark and quark-gluon splitting functions, for 
which the LO contributions does not include \/x terms and the resummation starts at NLL level. 
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Figure 3: The all-x gluon coefficient functions for the fragmentation functions Fj and Fi (multi- 



plied by x), down to extremely small values of the scaling variable x. As discussed above Eqs. (14.11) 
and (|5.1I) . the respective '(N)LO + resummed' approximations are defined by resumming as many 
logarithms as required to remove all l/x terms due to the (N)LO contributions. 
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Figure 4: As Fig. [3l but for the quark coefficient functions which are suppressed by one power of 
\nx w.r.t. the gluon quantities. As the diagonal splitting functions, Cj,q includes a 8(1 — x) term. 
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8 Summary 



We have presented the analytic all-order expressions in Mellin-A^ space for the next-to-next-to- 
leading logarithmic (NNLL) small-jc resummation of the splitting functions for parton fragmen- 
tation and the coefficient functions for the fragmentation functions Ft, Ft and F§ in gauge-boson 
and (heavy-top) Higgs -exchange semi-inclusive e + e~ annihilation. The resummation replaces the 
double-logarithmic x~ l \n n x small-* enhancement of the fixed-order results for all these quantities 
by an oscillatory behaviour which can be described in terms of Bessel functions for the splitting 
functions PZ (x, 0C S ) . The present results are sufficient to construct a combined 'NLO + resummed' 
approximations which is applicable down to extremely small values of x and determine the first 
three terms in the expansion of the first moments related to particle multiplicities in powers of y/Ok. 




In view of the dependence of the oscillation amplitudes on the logarithmic order and the large 
size of the N=l expansion parameter, an extension of the present results to the fifth (N 4 L) log- 
arithms and a 'NNLO + resummed' approximation is desirable. Completing these results is well 
beyond our present approach based on the D-dimensional structure of the unfactorized fragmen- 
tation functions and the mass -factorization relations. However, we have taken a first step in this 
direction by confirming the all-order relation between initial-state and final-state splitting functions 
suggested in Refs. ll2Tll22ll at the NNLL level and applying it to predict the 'non-singlet' Cp = 
part of Pj g at N 4 LL accuracy. 

A Form file of our NNLL results presented in the Sections 3, 4 and 5 above and Eqs. (IB.7I) - 
(IB. 101) below can be obtained by downloading the source of this article from the arXiv servers. 
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Appendix A 



The crucial step towards obtaining the analytic expressions presented above and in Appendix B 
was the solution of the sequence in Ref. |[T9l . The key to understanding this sequence was 
the observation that the denominators began as, and subsequently were divisors of, the sequence 
of least common multiples of triangular numbers [|36l . This suggested expanding each term as a 
sum of reciprocals of triangular numbers, which after some 'playing' with the initial terms led to 

A (2) _ J _ I 
qi 1 ~ 1 ' 

, (3) 112,3 

A qi ~ 3 _ 3 + 1 ' 

A (4) _ 73_5 7 9 

A qi ~ 6 _ 6 + 3 + 1 ' 



, (5) _ 1207 _ 14 , 19 , 23 , 28 
A qi ~ 30 ~ 10 + 6 + 3 + 1 



A (6) _ 2015 _42,56,66,76,90 rA n 

A qi ~ 15 _ 15 + 10 + 6 + 3 + 1 ' ••• {AA) 



and the recognition that the successive numerators are the coefficients of the series DTI which are 

(n) 

some particular sums of products of Catalan numbers. Our knowledge of A y qi ' to n = 17 checks this 
identification to the 136 th entries in the sequences [|36ll and 11371 . virtually excluding an agreement 
by coincidence. Using Maple, this result was rewritten in the more accessible form 



W _ 2(2n-2)! / 1 1 6 ,\ 2(2n)! 2 "y 3 1 
A qi ~ (n-l)!(n+l)! \n-\ + n + «+l £ ) + n!(n+l)! L k ' 



which was then understood to represent the expansion coefficients [in terms of ^ in Eq. (13.11) 1 of the 
first line in Eq. (13.21) . Once this function (the simplicity of which further supports the correctness 
of the above deductions) was known, it was not too difficult to derive all of Eqs. (13.21) - (13.51) . 

The corresponding x-space series can be readily obtained using Eqs. (IA.2I) and (17.11) . We have 
not been able to express the result as a known special function. It may be interesting for future 
studies, however, that the function L in Eq. (|3.1I) has a simple Mellin (and Laplace) inverse, viz 



[ dxxN ~ 2 \nx (M2^x) - l) = In (I + ^1 + ^ ) (A.3) 

with, in our case, a = SCacis- In view of Eq. (17.21) and the logarithmic jc-space arguments, the 
Mellin inverse of the first line of Eq. (|3.2I) and other products of the form (N — \ ) n (S z — 1)L with 
z= 1, —1, —3, . . . , can hence be expressed as Laplace convolutions of Bessel functions. 



Appendix B 

Finally we present various, mostly lengthy results beyond the (at N ^ 1) main 'NLO + resummed' 
approximation. We start with the N 3 LL corrections to the splitting functions Pi. and given at 
NNLL accuracy in Eqs. dill) and d3~3l) . 
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P q r q,N3L L W = 2552 ^f a s { ~ 3( 121 C A + MC\n f ~ ^C F n f + AClnj - \6C A C F nj 

+ \6C}nj) ^ ((S~ 2 — 1—4^— 16^ 2 ) - (S~ 3 - 1 - 6$)L) + 6(341^ + 400^ + SC%C F n f 
-AClnj + 32C A C F nj - A%cjnj) ^ (S^ 1 - 1 - 2\ - 6^ 2 ) + 2( [905 + 4608 X^c\ + 3456CjC F 
-\516Cln f + 86800^^ + 556C^n} - 2336C A C F /^ + 2064Cjn/) ^ (S - 1 + 2^ + 2^ 2 ) 
+3( [709 - 576 Xj\C\ - 300C%n f + 936C 2 C F n / - 16C%n} + A00C A C F nj - A96Cjnj) ^ (S" 1 
-1 - 2^)L - 24(55C^ - \2C\n f + 46C i 2 C F n / - 4C^n/ + 2QC A C F nj - 2ACjnj) ^ 
-1 - 2^)£ 2 - 6 ([5231 - 6048C 2 ]C A - [5184 - 6912^ 2 ]Q 3 C F - 896^^ - 2488C 2 C F n / 
-292Clnj + 2ll6C A C F n} - 3AAQCjnj) ^(S-1+2^)L- 24( [961 - 864 £ 2 ]C2 + [1 152^ 2 

-612]CjC F - 336C\n f + 3860^/^ + 4C 2 n / 2 + 92C A C F nj - 232Cjnj) ^ (5 - 1 + 2^)X 2 
-64( [6 + 18 £ 2 ]C^ + 21C\C F - 29C\n f + lOOC^C F n f + 8C>^ - 44C A C F n/ + 6QCjnj) 
• ^ (5 - l)L 3 + 192([175 - 288 £ 2 ]C^ - [261 - 288 fc] CjCp + A\C\n f -318C A z C F n f 

-2%Clnj + \16C A C F nj - 256Cjnj)^L 2 + 8([3659 - 9216 - [5832 - 8640 ^C^Cp 

-628Cjn / -3488C 2 C F n / -380C 2 ^ + 2560C A C F ^-3696C^)^(5- 1+2^) 



+ 12( [2017 + 576 + 1728CjC F - 1772Cjn f + 5A32C A 2 C F n f + 356C^nj 
-\lAAC A C F nj + 2Q6ACjnj) | (5 - 1 )£ - 384( [1 17 - 252 ^ 2 }C A - [207 - 288 £>] CjC F 

-\^ln f -l64C^C F n f -l4C}n} + 96C A C F n}-l44C}n})(l + ^L)^ (B.l) 

and 

<g,N3 LL W = §< q , N 3 LL (^) + i^ 2 {2([81-144yC A 4 -90CiC F + 144^ 2 CiC F 
-79Cjn f + \Q6ClC F n f + 6Cf^ - 2AC A C F n] + 24C*|n/) (5 - 1 + 2^ + 2^ 2 ) 
-2( [483 - 576 ^ 2 ]C A - [360 - 576 C 2 ]CjC F - 139C A 3 n f + 1 02(^0^ + 56C A C F nj 
-6Clnj - %%C}nj) | (S - 1 + 2^) + ( [429 - 576 ^ 2 ]C^ - [360 - 576 £ 2 ] CjC F - 2YhC\n f 
+250C A l C F n f - 2Clnj + 40C A C F n/ - 12Cjnj) ^ (S - 1 + 2£)£ + 8( [137 - 144 £ 2 ]C^ 
_ [90 - 144 C 2 ]C|C F - 17Q 3 n f - 6C^C F n f - 6C^nj + 36C A C F n} - AZCjnj) (1 + ^L) 
+(11C A * + \3C\n f - 26C A 2 C F n / + 2C|^ - 8C A C F n/ + ZCjn})^ (S" 1 - 1 - 2^)L 
+4( [59 - 72 £ 2 ]c£ - [45 - 72 ^ 2 ]CjC F - 20Cj{n f + 20C^C F n f - 2C A l n 1 f + \AC A C F nj 
-20Cjnj)^{S-\)L 2 Y (B.2) 

Together with the first moments of the corresponding truncated second-order splitting functions 
?q q ( V=l) = -(13-12^ 2 + 8^)C A C F + (26-24^ 2 + 16^ 3 )C^ + ^C F n / , (B.3) 
^qg ( V=l) = {^--^^C A n f + AC F n f , (B.4) 
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which can be derived, e.g., using Form packages [24J based on Refs. 113811391 as discussed in 
Ref. [19J, Eqs. (|B.1I) and (IB.2I) lead to the combined results 

^=1)^ = g {- [2106 - 1944 C 2 + 1296 C 3 ]Q 5 + [4477 - 2016 C 2 ]C A \ 

+ [4212 - 3888 £ 2 + 2592 £ 3 ] C\C F - \\2AC\nj + 1728 & CjQ^ 

+ 116C2nj+2336CiCFn/-448C A C F n/+432C|nj| , (B.5) 

P^(N=l)\ a2 = ig2^ {[3913 - 1 152 £ 2 ]Cl+ 1728 CjC F - 812^+ \16Clnj 

+22A0C A 2 C F n f - 5AAC A C F nf + 624 C|n^} . (B.6) 

The numerical values for ny = 5 of these coefficients have been included in Eq. (13.81) above. 

The NNLL corrections to the coefficient functions (14.11) - (14.41) for Ft and F§ can be written as 

Cr, q ,NN L (A0 = 452^?«s{-192(5C A 2 -4C A n / + 12C F n / )(^i: + l)+8(35C A 2 -14C A n / 
+24C F n / )|(F~ 3 -1 + 35) -8(35C A 2 + 34C A n / - V2QC F n f )\{F- 1 - 1 +5) -3 (505 Cj 

- 150C A n / + 416C F n / )|(F" 1 - 1)X — 4(887 C A — 122C A «y + 336C F ny)(F — 1) 

+ 4 ( 1 69 + 86 Caw/ - 276 C F n f ) | (F - 1 ) £ + 1 2 ( 1 1 C| + 2 Qny - 4 C F /^) (4 (F 3 - 1 ) 

-i(F 5 -l)iL) + (llC A 2 -2C A n / )(20(F 5 -l)+6|(F 3 -l)i:-5i(F 7 -l)i:) 

+ 96(5C A 2 -2C A /i / + 6C F n / )^Fi: 2 | (B.7) 

CtannlM = -^C rq;NNL (^V) + g2544^a s {9216([54-72^ 2 ]C A 4 - [45 -72^]C A 3 C F 
-23 Cjn 7 + lOClCprif) (1 - 1 (5 - 1 + 25 - X)) - 36864 (C%C F n f + C A C F nj 
-2 Cjnj) I ( (5- 1 + 2^) + (F- 1)X) +256 ([637- 1620 [810- 1296 ^ 2 ]C A 3 C F 

- 206Cln f + 129(^0^ - 44Cjv^ + 258 C a C f h 2 - 342 C|n/) ^ (F~ 3 — 1+3^) 
+ 256 ([1523 - 972 £ 2 ]c£ - [810 - 1296 ^ 2 ]C|C F - 2<&lC\n f + 159(^0^ - 28 

+ 318 C A C F n 2 - 522 C 2 F nj) ± (F~ 1 - 1 +£) + ([1285825 - 1327104^ 2 ]C A - 331776 
■ [1 - 2 C 2 ]CjC F - 2Al\6Cln f - 19392C A 2 C F n f - 45020C A 2 nj + 233088C A C F nj? 

- 290304 C}n 2 f ){F - 1) + 8 (14839(7^ + 1 168 Cj/^ - 2\12ClC F n f - 2%AC\n 2 f 
+ n2SC A C F nj -\\52Cjnj){F 3 - 1) + ([75461 + 165888 £2]^ + 41 1560^ 

- 1 \00\6ClC F n f + 1 UAAClnj - 5A\AAC A C F nj + 6220%C}n}) (F 5 - 1) + 16 (1 1957 C\ 

- A0Cln f - 4488C A 2 C F n / - 388C^ + 816C A C F n / 2 )(F 7 - 1) - (21901 C\ + 55396^^ 

- 63360C^C F n f + 724Cj^ - 1 \520C A C F nj + 1 1520C jnj){F 9 - 1) - 840 (121 c\ 

- 44 ClC F n f -4C%n} + 8 C A C F nj) (F 1 1 - 1 ) - 3 85 C\ ( 1 1 C A - 2 n f ) 2 (F 13 - 1 ) j (B .8) 
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q,, g ,NN L (AO = -C rq;NNL (AT) + ^44 {256 ([601 -32ACAC} + 64C}n f + 93C}C F n f 
-%Clnj - 30C A C F nj + 90Cjnj) ^ ((S - l)F~ l + 2$) - 1 8432 (2C%C F n f - C A C F nj 
+ 2C}nj) ^(S-l)FL+ ([3335233 - 663552 £>]Ca - 162380^^ - 356928C A 2 C F n f 

- IQQAClnj + A6A6AC A C F n] - 69\2QCjnj) (F - 1) + (9559 C\ + 3256(7^ + 244 c\n} 

- 6732C%C F n f - \Q%QC A C F nj + \\52C}nj) 8 (F 3 - 1) - ([55483 - 165888 Ca]C} 
-175A0C A 3 n f + 52032C A 1 C F n f -AlA&C A l n} + 30336C A C F n}-A3776C F l n})(F 5 -l) 

+ 16 (9317C£ - (220C A 3 + 3l68C A l C F )n f - (268C^ - 516C A C F )rf)(F 7 - 1) - (21901C^ 
+ 55396Cj«y - 63360C}C F n f + 12AClnj - 1 \520C F nj(C A - C F )) (F 9 - 1) - (121 C A A 
-AAClC F n f -ACln] + %C A C F n})%AO(F n - 1) - 385 (1 IC A - 2n f ) 2 Cj((F 13 - 1) j (B.9) 

q,, q ,NN L (A0 = ^C Tq)NNL (N) + ^a s {-2(C A 2 + C A n f -2C F n f )(l + ^L) 

- (2C A 2 -C A n f + 2C F n f )± ((F~ 3 - 1 + 3$) - (F" 1 - \)L) 

+ (AC 2 + C A n f -2C F n f )^ ((F _1 - 1 + £) - (F - 1)X) } (B.10) 

Finally we write down the N 4 LL contribution to the splitting function P^ g in the limit C F = 0, 
which still contains one unknown coefficient as discussed above, 

^WISll = 1327 1 1040 ^ a s 2 ^{l6([15688235- 19918080^ + 7983360^3 

+ 5059584 C|]Ca + [914360 + 875520 £ 2 - 2142720 ^Cj^ -[134200 + 46080 ^ 2 ]C A n} 
+ 5600njc A - 80n/) (5 - 1) - 32( [7822505 - 2826000 & + 1330560 £ 3 - 134784 £f ] c\ 
+ [5 14490 + 83520 ^ 2 - 276480 £ 3 ] Cj W/ + [16880 + 20160 £ 2 ] C^n/ - 2840 wjC A + 80 n/) 

• | (5- 1 +2£) + 2( [12686895 + 12997440 £ 2 - 247 1040 £ 3 - 10907136 X£\c\ 

- [3309880 + 564480 £ 2 - 1624320 £ 3 ] C\n f + [37960 - 172800 &}C A nj + 21280 n 3 f C A 

- 1040/iyf) (5 - 1 + 2^ + 2^ 2 ) - 4( [3 135445 + 6822720 £ 2 + 1 90080 £ 3 - 5868288 £f ] 

- [1973 120 + 587520 ^ 2 - 1071360 £s] C\n f + [106520 - 149760 £ 2 ] C^n/ + 14080njC A 

- 1200n/) ^ (S^ 1 - 1 - 2^ - 6^ 2 ) - ( [2095591 + 158976 & - 1 140480 £ 3 + 33 1776 £f ] 
+ [61560 + 396288 - 207360 £ 3 ] C A «/ - [83352 - 64512 ^ 2 ] C^n/ + 224 n 3 f C A + 880 n/) 

• 5 ^ (S~ 3 - 1 - 6^ - 30^ 2 ) - 10( [198803 - 209088 £ 2 ] C A + [69872 - 76032 £ 2 ] Cj^ 

- [600 + 6912 y Clnj - 236%n 3 f C A - 208 n}) h (S~ 5 - 1 - 10^ - 70^ 2 ) 

-25(llC A + 2n / ) 4 ^(5- 7 -l-14^-126^ 2 )} + Cja s 2 /Y (S + S' 1 -2)BgP , (B.ll) 
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where the normalization of the last coefficient has been given below Eq. (16.8b . This result leads to 
the first moment 

Pg g (N=l)\ C a f° = (2C A ) 1 / 2 4 T ^ 7 2^{[182872 + 175104^ 2 -428544^]C A 3 n / 

+ [3137647-3983616^2+ 1596672^ + 5059584/5 Q]C% (B.12) 
- [26840 + 9216 £ 2 ] C%nj+ U20C A rf - 16nf + 165888C|fl2 3) } 

which has been included in an approximate numerical form for n^ = 5 Eq. (16.81) . 
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